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Abstract 

We propose that the radiative decay process, % — * J/ip 7, may serve as a clean searching mode for 
ijb in hadron collision facilities. By a perturbative QCD calculation, we estimate the corresponding 
branching ratio to be of order 10~ 7 . Though very suppressed, this radiative decay channel in fact 
has larger branching ratio than the hadronic decay process % — > J/tpJ/ip, which was previously 
hoped to be a viable mode to search for r]b in Tevatron Run 2. The discovery potential of rjb through 
this channel seems promising in the forthcoming LHC experiments, and maybe even in Tevatron 
Run 2, thanks to the huge statistics of % to be accumulated in these experiments. The same 
calculational scheme is also used to estimate the branching ratios for the processes % (r/ c ) — * (fry. 
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I. INTRODUCTION 



The existence of rjb, the pseudo-scalar partner of T(1S), is a firm prediction of QCD, about 
which nobody would seriously challenge. It is rather unsatisfactory that although three 
decades have elapsed since the discovery of T(IS'), this particle eludes intensive experimental 
endeavors and still eagerly awaits to be established. 

On the theoretical side, many work have been devoted to uncovering various properties 
of r]b, such as its mass, inclusive hadronic and electromagnetic widths, transition rates and 
production cross sections in different collider programs Among all the observables, its 
mass is believed to be the simplest and most unambiguous to predict. Recent estimates for 
T — 7/6 mass splitting span the 40-60 MeV range 0,13, Q, [5]. An eventual definite sighting 
of rjb and precise measurement of its mass will critically differentiate varieties of theoretical 
approaches, consequently sharpening our understanding towards the bb ground state. 

The 776 has recently been sought from 77 collisions in the full LEP 2 data sample, where 
approximately few hundreds of rjt, are expected to be produced. ALEPH has one candidate 
event with the reconstructed mass of 9.30 ± 0.03 GeV, but consistent to be a background 
event 6]. ALEPH, L3, DELPHI have also set upper limits on the branching fractions for 



7/6 decays into 4, 6, 8 charged particles 



8jl. Based on the 2.4 fb 1 data taken at the 



T(25) and T(3S) resonances, CLEO has searched distinctive single photons from hindered 
Ml transitions T(2S),T(3S) — > rj^j, and from the cascade decay T(3S) — > h b n°, h b n + n^ 
followed by El transition — > 7/67, but no signals have been found [9j]. 

Hadron collider experiments provide an alternative means to search for 7/6. Unlike the 
e + e~ machines which are limited by the low yield of 77^, hadron colliders generally possess a 
much larger 7/6 production rate, which in turn allows for triggering it through some relatively 
rare decay modes yet with clean signature. However, one caveat is also worth being em- 
phasized. The corresponding background events are in general copious in hadron machines, 
so the virtue of this sort of decay modes may be seriously discounted (Such an example is 
the electromagnetic decay 776 — > 77, with an expected branching fraction ~ 10 -4 , which is 
nevertheless overshadowed by the ubiquitous 7 events originating from 7r° decay). 

Several years ago, Braaten, Fleming and Leibovich suggested that the hadronic decay 
Vb J/ipJ/ip, followed by both J/ip decays to muon pairs, can be used as a very clean 
trigger to search for 776 at Tevatron Run 2 lOj. By some simple scaling assumption, they 
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estimate the branching ratio of the double J/ip mode to be 7 x 1CT 4±1 , and conclude that 
the prospect of observing the r\\, — > 4/j channel at Tevatron Run 2 is rather promising. CDF 
in fact has followed this suggestion and looked at the full Run 1 data for the 4/x events in 
the expected r)b mass window 11]. 

However, some objection has been recently put forward by Maltoni and Polosa, who argue 
that the estimate of the branching ratio of double J/ip mode by Braaten et al might be too 



optimistic 



121 ] . They instead advocate that rjb — » D*D , with the estimated decay ratios 



lying in the 10~ 3 — 10~ 2 range, may serve as better searching modes for 775 in Run 2. 

Very recently, one of us (Y. J.) has surveyed the discovery potential of various hadronic 
decay channels of r\ h jis] ]. The explicit perturbative QCD (pQCD) calculation predicts 
Br[r/5 — > J/ipJ/ip] to be only of order 10~ 8 . If this is the case, the chance of observing this 
decay mode in Run 2 then becomes rather gloomy, whereas the observation possibility at 
LHC may still remain. Another noteworthy assertion of 13| is that, by some rough but 
physical considerations, one expects Brfat -> D*D] ~ 1(T 5 and Br]^ -> D*D } ~ 1(T 8 , 
which are much smaller than the optimistic estimates by Maltoni and Polosa. Taking the 
low reconstruction efficiency of D mesons further into consideration, these double charmful 
decay modes may not be as attractive as naively thought. 

In this paper, we try to propose another decay process, T]b — > J / ^7, as a viable discovery 
channel for rjb in hadron collider experiments. At first sight, this mode, being a radiative 
decay process, may not look very economical due to the large total width of rjb- In fact, 
our explicit pQCD calculation reveals that the corresponding branching ratio is indeed very 
suppressed, only about 10~ 7 . Nevertheless, it is worth emphasizing that this is already larger 
than that of the double J/ip mode. The usual way of reconstructing J/ip is via the dimuon 
mode, with a branching fraction about 6%. In the light of this, the advantage of this radiative 
decay process over the double J/ip mode will be further amplified, since only one J/ip needs 
to be reconstructed in our case. Unlike the open-charm decay processes r\b — > D*D , in 
which the final decay products such as K, 7r in general suffer severe contamination from the 
combinatorial background, the presence of J / ip in the final state renders our radiative decay 
channel much cleaner to look for experimentally. 

Thanks to the huge amount of rjb to be accumulated in high energy hadron collision 
facilities, we expect that this radiative decay channel, albeit being a rare decay mode, may 
still have bright prospect to be observed in Run 2 and in the forthcoming LHC program. 
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FIG. 1: Lowest-order diagrams that contribute to % — > J/ipj. Diagrams (1), (2), (3) stand 
for QCD-initiated process, while (4) represents the dominant QED process governed by photon 
fragmentation. Crossed diagrams are implicitly implied. 



Not surprisingly, one should be aware that the major QCD background events, i.e. the 
associated J/ip + 7 production [ij], 16], may exceedingly outnumber our signal events. 



To make this mode practically viable, one has to ensure that those background events can 
be significantly depressed by judiciously adjusting the kinematical cuts. 

The remainder of the paper is distributed as follows. In Sec. flU we present the pQCD 
calculation for the radiative decay process rjb — > J/ipJ, treating heavy quarkonium states 
in NRQCD approach which, to our purpose, is equivalent to the color-singlet model. In 
Sec. IHIt we present the numerical prediction to the corresponding branching ratio, and 
explore the observation potential of this decay mode in Tevatron Run 2 and in the coming 
LHC experiment. We subsequently apply the same formalism to estimate the branching 
ratios for analogous processes rjb (t) c ) ~^ 4>l-> which is equivalent to work in the context of 
constituent quark model. We summarize in Sec. HVl 



II. COLOR-SINGLET MODEL CALCULATION 



In this section, we present a pQCD calculation for the decay rate of 77?, — > J/iftj. This 
decay process can be initiated by either strong or electromagnetic interaction, with the 
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corresponding lowest-order diagrams shown in Fig. [TJ Since the annihilation of the bb pair, 
the creation of the cc pair, as well as the emission of the hard photon, are all dictated by 
short-distance physics, with the hard scales set by the heavy quark masses, it is appropriate 
to utilize the pQCD scheme to tackle this exclusive process. 

Nowadays it has become standard to employ the NRQCD factorization framework to 



cope with hard processes involving heavy quarkonia [17J. One of the major advantages 
of this model-independent EFT approach over the old color-singlet model is that it can 
in principle incorporate the contribution of higher Fock components of quarkonium (color- 
octet effect) in a systematic fashion. This is exemplified by the systematic treatment of 



inclusive annihilation decay of various quarkonium states within this framework [171 . Il8| . In 
contrast, the inclusion of color-octet effect in exclusive processes has not been developed to 
a comparable level within this context. To date, there are only few cases where this effect 
has been investigated in the model-independent language. One example is the magnetic 
dipole transition process such as J/ip — > rj c ^, where the leading color-octet contribution 



vanishes |lfl| . Another example is tin* quarkonium radiative decay to light mesons, such as 



T — > /2(1270)7, where the color-octet effect turns out to be quite insignificant 20]. 

If we take the lessons from the aforementioned radiative decay processes, it seems per- 
suasive to assume that the color-octet effect may also be unimportant in our case. We will 
completely ignore this effect in this work. In this regard, there will be no difference between 
the NRQCD approach and the color-singlet model, consequently these two terms will be 
used interchangeably. 

Before launching into the actual calculation, it is worth mentioning that, Guberina and 
Kuhn has studied the analogous radiative decay process T — > 7^7 in NRQCD approach more 



than two decades ago 211 ] . Our calculation will intimately resemble theirs. 

It is useful to take notice of a simple trait of this process, that the photon can only be 
transversely polarized, so is the recoiling J/ip by angular momentum conservation. This 
property holds true irrespective of whether this process is initiated by strong or electromag- 
netic interaction. In fact, parity and Lorentz invariance constrains the decay amplitude to 
have the following unique tensor structure: 

M{\i, A 2 ) = Ae fWaP £f /ij {X 1 )E; v {X 2 )Q a k . (1) 

We use Q, P and k to signify the momenta of rj b , J/ip and 7, respectively, and use Ai 



and A2 to label the helicities of J/ip and 7 which are viewed in the rjj, rest frame. It is 
straightforward to infer from (JT]), that the only physically allowed helicity configurations are 
(Ai, A2) = (±l, ±1)- All the dynamics is encoded in the coefficient A, which we call reduced 
amplitude. Our task then is to find out its explicit form. 

In the color-singlet model, it is customary to start with the parton process 6(p&) b(pi) — > 
c(p c ) c(pc) + then project this matrix element onto the corresponding color-singlet 

quarkonium Fock states. For reactions involving heavy quarkonium, it is conventional to 
organize the amplitude in powers of the relative momentum between its constituents, to 
accommodate the NRQCD ansatz. This work is intended for the leading order accuracy 
in relativistic expansion only, it is then legitimate to neglect the relative momentum inside 
both rjb and J/ijj. We thus set pb = p^ = Q/2 and p c = Pc = P/2. For the bb pair to form 
r] b , it is necessarily in a spin-singlet and color-singlet state, and one can replace the product 
of the Dirac and color spinors for b and b in the initial state with the projector 

u( Pb )v(p- b ) — > ^=($ + 2m 6 )* 75 x ^J=^(0)J ® (-^=J . (2) 

For the outgoing J/ip, one can employ the following projection operator: 

v(p,)u(Pc) — ^= f J/il {P + 2m c ) x ^_^^(0)J ® (-^=J , (3) 

where £j,^ is the polarization vector of J/ip satisfying Ej/^(\)-e*j^(\') = —5 XX ' and P-e = 0. 
A^ c = 3, and l c stands for the unit color matrix. The nonperturbative parameters, Vv(0) 
and ipj/^(0), are Schrodinger wave functions at the origin for rjb and J/ip, which can be either 
inferred from phenomenological potential models or directly extracted from experiments. By 
writing (j2J) and ([3]) the way as they are, it is understood that M Vb = Im^ and Mj/^p = 2m c 
have been assumed. 

We commence with the strong decay amplitude. As can be seen in Fig. (U the lowest 
order QCD-initiated contribution starts already at one loop. Note that the photon can only 
be emitted form the c quark line, because the even C-parity of rjb forbids the 7 to attach to 
its constituents. Using the projection operators in (T5]) and ([3]), it is straightforward to write 
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down the Feynman rules for the strong decay amplitude: 
-Mstr = -C str e c eg 



4^(0)^(0) fd A h 1 1 



8^mWT c J (2vr) 4 k\ k\ 
trjjjg + 2m b )^u{A- h + m b )j Pi , , , , 

X — — 2 r yjJ' *^ V, K\ <— > K 



(Pb ~ h) 2 - m 



2j 



X 



+ 



tr[^(f + 2m c ) 7 ^ c - ft + m c ) - j4 + m c )Y\ 



{{p c - h) 2 ~ ml) ((p e - /c 2 ) 2 - m 2 c ) 
tr[fj^(P + 2m c ) f^ c + }6 + m c ) 7 ^ 2 - f- c + m c ) 7 "] 



((p c + - ml) ((p s - A^) 2 - m 



2^ 



tr[^( > P + 2m c )7^(^ c - fc + m c )7"(- ^- # + m c ) 

+ ((pc-hy-mDdpt + ky-mi) < { ] 

where the corresponding color factor C stT = iV" 1 tr(T a T b )tr(T a T b ) = |. The momenta 
carried by two internal gluons are labeled by ki, A; 2 , respectively, which are subject to the 
constraint k\ + k 2 = Q. Notice that the box diagrams in Fig. [1] are related to one-loop four- 
or five-point functions, and the corresponding loop integrals are ultraviolet finite. Moreover, 
the occurrence of heavy b and c masses ameliorate the infrared behavior of the loop integral 
so that the result turns out to be simultaneously infrared finite. Since there is no need for 
regularization, we have directly put the spacetime dimension to four. 

After completing the Dirac trace in (Tj0), we end up with terms which do not immediately 
possess the desired tensor structure of ([I]), instead with one index of Levi-Civita tensor 
contracted to the loop momentum variable. Of course, when everything is finally worked out, 
all these terms must conspire to arrive at the desired Lorentz structure. Conversely, one may 



exploit this knowledge to simplify the algebra prior to performing the loop integral 2ll. 122]. 
First pull out the partial amplitude through M. s tr = e*f^ e* u . Eq. (JTJ then demands 

M„ u = Atr e^Q a k^ (5) 

which is compatible with the expectation that would vanish unless fi, v are in trans- 
verse directions. Contracting both sides of (jSJ) with e^^Qpk^, we can express the reduced 
amplitude as 

At r = ^^y 2 M, v Q p K . (6) 

After this manipulation is done, it is convenient to adopt a new loop momentum variable 
q, which is related to the old one via k\ = (Q + q)/2 and k% — (Q — q)/2. We end in a 



concise expression: 



6 7T 2 



m c ^(0)^/^(0) /^m 




m;, (mf — ml) 2 



mi 



(?) 



where 



(fc-Q) q 2 -{k-q) (Q-q) 



(8) 



»7T- ,/ - (q + Q) 2 (q - Q) 2 (q 2 + 2k-q- P 2 ) (q 2 - 2k ■ q - P 2 ) 

Since / is dimensionless, it can depend upon m& and m c only through their dimensionless 
combination. It is interesting to note that the b quark propagator has been canceled in this 
expression. The loop integral can be performed analytically by the standard method, and 
the result is 
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(9) 
(10) 



2 — it ' 1 + M 

where Li 2 denotes the dilogarithm (Spence) function. We will expound the derivation of this 
result in the Appendix. Note / has an absorptive part, which reflects that the intermediate 
gluons are kinematically permissible to stay on shell. One can apply Cutkosky's cutting 
rule to verify (ITUI) . The real and imaginary parts of / as function of m 2 c jm\ are displayed 
in Fig. 2. 

It is instructive to know the asymptotic behavior of / in the u — > limit, which can be 
readily read out from Qj and (ITU]) : 



Ref(u) 
Im/(tt) 



7l~ 



(1-2 In 2) lnw+ — -ln2 + ln 2 2 
3 

?r(l- 2 In 2). 



(11) 
(12) 



Note as u approaches 0, the imaginary part remains finite, whereas the real part blows up 
logarithmically. This trend can be clearly visualized in Fig. 2. The logarithmic divergence 
in the m c — > limit is obviously of infrared origin. Nevertheless, this does not pose any 
practical problem, since a non-relativistic description for a zero-mass bound state, as well 
as the resulting predictions, should not be trusted anyway. It is interesting to note that, the 
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FIG. 2: Real and imaginary parts of f(u) 



analogous radiative decay process, T — > 7^7, has a different asymptotic behavior, in which 



2l| 



both the real and imaginary parts of the QCD amplitude admits a finite limit 

We next turn to the pure QED contribution to this radiative decay process. One may 
naively expect that this contribution is much less important than the QCD contribution. 
However, it turns out that this expectation is not true and the QED contribution must be 
retained. Due to the neutral color charge of photon, this radiative decay can arise at tree 
level from the photon fragmentation, as shown in Fig. [1] (4). Obviously, the fragmentation 
type contribution is much more dominant over other types of QED diagrams. To calculate 
this contribution, a necessary input is the 7 — J /if) coupling, which is characterized by the 
J /if) decay constant 1 : 

(J/if)(X)\^c\0) = - gj/i) e^(X). (13) 

In the non-relativistic limit, gj/^ is linked to i/jj/^(Q) through the relation gj/^ = 
2^ 2 Nl l2 m l J 2 i)j/^), which can be derived from ([3]). The calculation is much easier than 



1 In conformity with the sign convention adopted in ([3]), a minus sign is compulsory to put in here, which 
takes into account the Grassmann nature of the quark field operator. 
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its QCD counterpart, and the reduced QED amplitude is 

^ ^(0)^/^(0) 



A 



e 3 iV n 



771* ) 



Substituting (GO) and ( 1141) into the formula 



47T 



| Atr + ^ 



(14) 



(15) 



we then obtain the desired partial width. Here |k| = (m| — m 2 c )/m b is the photon momentum 
in the r] b rest frame. This formula already takes into account the sum over transverse 
polarizations of both J /if) and 7. 

For phenomenological purpose, it is instead more convenient to have an expression for 
the branching ratio, where if) Vb (0) drops out: 



Br [7ft -> J/ipi] 



'dotal 



m r 



3/T ml (ml — ml) 



nn 



m; 

9 I — I 

mi 



where 



9(u) 



9 ii el a 1 



it 



2ai 



u 



(16) 



(17) 



encodes the electromagnetic contribution. In deriving this, we have approximated the total 
width of r\ b by its gluonic width: 

Sua 2 



3 m£ 



<(0) 



(18) 



where the LO expression in a s and v b is used for simplicity. 

Eq. ffTB"]) constitutes the key formula of this work. This equation conveys that, despite 
the adversity caused by the suppression a /a 2 , the QED fragmentation contribution never- 
theless enjoys the kinematic enhancement of ml jm\ relative to the QCD amplitude. For 
the physical masses of b and c, these two competing effects have comparable magnitudes. 
In the asymptotic regime as m b /m c ^> 1, the QED amplitude will eventually dominate over 
its QCD counterpart, because the enhancement factor ml/ml of the former is much more 
eminent than the mild \og(ml/m 2 c ) rising of the latter. At any rate, it is imperative to in- 
clude the QED fragmentation contribution. Moreover, it is important to recognize that the 
interference between these two amplitudes is destructive. This is opposite to what occurs in 

Vcl 0] and in the decay r/ b — > J/ijj J /if) [3], where the 



the analogous decay process T 
interference is constructive. 
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III. PHENOMENOLOGY 



A. Observation Potential of % — > in Tevatron and LHC 

It is now the time to explore the phenomenological implication of ([TBI) . The input pa- 
rameters are mj, m c , a, a s and ipj/^(0), all of which can be inferred from other independent 
sources. The wave function at the origin for J/ ip can be extracted from its electronic width: 

Air p 2 rv 2 

r[j/^e+ e -] = -^-V^(o), (19) 

where the LO formula in a s and v 2 c is used for simplicity. Using the measured dielectron 
width 5.55 keV [23J, we obtain Vj/V>(°) = °- 205 GeV 3/2 for m c = 1.5 GeV. Taking m b = 
M Vb /2 « 4.7 GeV, m c = 1.5 GeV, a = 1/137 and a s (m b ) = 0.22, we then find 

Bi[r] b -> = (1.5 ±0.8) x 10~ 7 . (20) 

The uncertainty is estimated by varying a s (/i) between 0.18 and 0.26 (which corresponds 
to slide the scale from /x = 2m b to 2m c ), as well as taking into account the errors in 
the measured T e + e - (of ±0.14 keV). The destructive interference between electromagnetic 
and strong amplitudes has pronounced effect. For the central values of input parameters, 
omitting the QED contribution will result in a prediction of 3.5 x 10~ 7 , which is more than 
twice larger than the actual value. To develop a concrete perception, we enumerate the 
values of / and g evaluated at u = m c /ml = 0.10 and a s = 0.22: 

/ = 3.5 -2.4i, = 2.1. (21) 

Clearly, Re/, Imf and g all have comparable magnitudes. As a result, the destructive 
interference effect is particularly important. 

The numerical prediction presented in ( J20l) is obtained by only using the tree-level match- 
ing coefficients for the total r\ b width and leptonic width of J/ip. One may worry that this 
oversimplified procedure will induce some error because it is known that the next-to-leading 
perturbative corrections to both quantities, especially to the J/ip leptonic width, are large. 
Let us assess their effects now. To the NLO accuracy in a s , one needs to multiply Eq. (TlBl) by 
1 ± (53/2 — 3l7r 2 /24 — 8n//9) a s (2m b )/^ (nf stands for the number of active light flavors), 
as well as multiply Eq. (TT9T) by (1 — 8a s (2m c )/37r 

) 2 H- 

Incorporating these corrections 
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amounts to multiplying (I20j) by a factor 2 



10.2a s (2m b )\ 1 / 8a s (2m c NX 2 



1 + ^-^^J ~ 1-04, (22) 

where n/ = 4 has been taken. Since including the perturbative radiative corrections has 
negligible net impact, we will keep using (120]) in the following phenomenological analysis. 
It is enlightening to compare (|20|) with the NRQCD prediction to the branching ratio for 



rjb decays to double J/if> [l3| |: 



Bi[7] b -> J/ip J/ip] = 2.4j£g x 10~ 8 . (23) 

Notice that the branching ratio of our radiative decay process is almost one order-of- 
magnitude larger than that of this hadronic decay process! The reason can be traced as 
follows. The double J /if} decay mode, though being a hadronic one, has maximally violated 
the helicity selection rule of pQCD. As a result, the branching ratio gets severely suppressed, 
oc a 2 s vl (m c /mb) s [3]. In contrast, if we count / ~ 0(1) 3 , Eq. (fIBl) then implies that our 
radiative decay process admits a scaling behavior Br ~ a a 2 s v% (m c /mb) , which is much 
more mildly suppressed by powers of 1/m 2 . and v c relative to rjb —>■ J/ipJ/ip, hence is more 
favorable. 

Experimentally J /if} can be cleanly reconstructed by decays to lepton pairs. Multiplying 
( 120]) by the branching ratios of 12% for J /if) decays to and e + e~, we obtain Br[rjb — > 

J/ip'-f — > l + l~ r f\ = (0.8 — 2.8) x 10 -8 . The total cross section for producing r/b at Tevatron has 



been estimated to be about 2.5 /zb [12]. The production cross section for the events 
is thus about 0.02 — 0.07 pb. For the full Run 1 data of 100 pb _1 , we then obtain between 
2 and 7 produced events. Because the kinematical cuts, as well as taking into account the 
acceptance and efficiency for detecting leptons, will further cut down this number, it seems 
not so fruitful to assiduously seek the rj b through this radiative decay mode in Run 1 data 
sample. 

Tevatron Run 2 aims to achieve an integrated luminosity of 8.5 ftT 1 by 2009. Assuming 
equal a(pp — > r)b + X) at \fs = 1.96 and 1.8 TeV, we then estimate there are about 200-600 



2 It should be kept in mind that the relativistic correction to leptonic decay of J ftp is also large. We have 
not considered this complication. 

3 The slow rising \u{m\ / m 2 c ) term in Re/ is of no concern here for physical b and c masses. Also we neglect 
the pure QED contribution for the lucidity of the argument. 
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produced events. The product of acceptance and efficiency for detecting J/ if) decay to muon 
pair is estimated to be e « 0.1 [lOj. It may sound reasonable to assume the corresponding 
factor for the electron pair also of the same magnitude. Multiplying the number of the 
produced events by e, we expect between 20 and 60 observed events in the full Run 2 
period. This is quite encouraging, but we should be cautious about the fact that the major 
QCD background events, the associated J/ift + ^y production, could also be copious. At the 
HC colliders, the dominant production mechanism for these events is through 



Tevatron anc 
gg fusion [l4 



161 ]. A naive analysis indicates that the direct J/ if) + 7 production with 



an invariant mass near 9.4 GeV preponderates over that from the rjb decays. 

It is important to keep in mind that the J /if) stemming from the radiative 77^ decay 
must be transversely polarized. This characteristic may be used to discriminate the signal 
events from the background events, because J /if) from the latter processes can also be 
longitudinally polarized. Furthermore, the detailed kinematical distributions of background 
processes need to be thoroughly studied, in order to guide experimentalists to choose the 
optimal kinematical cuts to oppress as many background as possible, and in the meanwhile 
without significantly sacrificing the signal events. This is somewhat beyond the scope of this 
work and needs further independent studies. 

The forthcoming LHC experiments will greatly increase the number of the produced 
775 — > l + l~ r y events. To assess the discovery potential of this mode at LHC, we need first 
know the inclusive production rate for 775. There are rough estimates for the Xbo,2 cross 
sections at LHC, which are about 6 times larger than the corresponding ones for producing 
them at Tevatron j^. Assuming the same scaling also holds for rjb, we then obtain the cross 
section for r] b at LHC to be about 15 fib, subsequently the production cross section for the 
/ + /~7 events to be about 0.1-0.4 pb. For 300 fb~ data, which is expected to be collected 
in one year run at LHC design luminosity, the number of produced events may reach about 
3xl0 4 — lxlO 5 . Multiplying the number of the produced events by e, we expect between 
3 x 10 3 and 1 x 10 4 observed events per year. 

Based on this analysis, we are tempted to conclude that, the chance of observing rjb at 
LHC through this mode is very promising. With such a large amount of signal events, 
it is possible to measure the leptonic angular distribution to pin down the polarization of 
J /if). As has been stressed, in order to effectively selected the signal events out of the 
abundant background events, one needs to develop a thorough understanding towards the 
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QCD background. 



B. Radiative decay of rj b (ry c ) into 



When coping with light mesons in hard exclusive processes, the most natural treatment for 
them is using light-cone wave functions. Indeed, heavy quarkonium radiative decays to light 
mesons, such as T — > / 2 (1270)7 20|, |25] (Ref. |20|] has heavily exploited the EFT machinery, 



i.e. NRQCD combined with SCET), J/tp(Y) — > r]{r]')^ [26j, |27|, have been studied along 
this line. On the other hand, the constituent quark model, which treats the light meson as a 
non-relativistic bound state, is also frequently invoked as an alternative method for a quick 
order-of-magnitude estimate. Various radiative decay processes, e.g. J ftp decays into light 
pseudo-scalar and P-wave mesons, as well as Xcj into light vector mesons, have already been 



studied in this context 



22 



28]. 



We may apply the same strategy to analyze the radiative decay t] b {rj c ) -> K7. We take 
Vb ^ 07 as a representative. By regarding as a strangenium, we can directly use (|16p . only 
with some trivial changes of input parameters. We take the m s = « 0.5 GeV. The 

wave function at the origin of 0, t/u(0), can be extracted through (Tl9|) from its measured 



electronic width of 1.27 ± 0.04 keV 
constant between a s {2m b ) = 0.18 and a s {2m 
uncertainty in T e + e -, we obtain 



231 ] . Taking mj = 4.7 GeV, varying the strong coupling 
= 0.51, and including the experimental 



Br[r] b -> 7 ] = (0.3 - 6.9) x 10 



-8 



(24) 



The QED contribution dominates in this case due to the larger ratio of m b to m s . To see 
this concretely, we list the values of / and g evaluated at u — m 2 s /ml w 0.01 and a s = 0.22: 



./ 



4.9 - 1.5i 



g = 23.5. 



(25) 



The dominance of g over |/| is apparent. Neglecting QED contribution for this a s value 
will decrease the branching fraction by one order of magnitude. For the absence of clean 
signature for 0, such a rare decay mode will be extremely difficult to observe in hadron 
collision experiments such as LHC It may be interesting to compare this radiative decay 
channel with the following hadronic mode 131 ]: 



Br[r] b 



(0.9-1.4) x 10" 



(26) 
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This 20 mode is even smaller, and hopeless to be seen experimentally. 

One can proceed to consider the analogous decay process r\ c — > 07. Parallel to the 
preceding analysis, taking m c = 1.5 GeV, and varying a s (/i) from a s {2m c ) = 0.26 to 
a s {2m s ) = 0.51, we obtain 

Bi[r] c -> 7 ] = (2.1 - 8.6) x 10~ 7 . (27) 

The destructive interference pattern is similar to that in rji, — > J/V'T- This decay mode is 
still too much suppressed to be observed in the functioning charmonium factory like BES II, 
perhaps also in the forthcoming BES III experiment. It seems also rather difficult to observe 
this decay mode in the current and future hadron collision facilities. 

IV. SUMMARY 

The motif of this work is to suggest one viable way to ferret out i] b in the functioning 
and forthcoming hadron collider facilities, that is, through its radiative decay into J/ip. 
This decay mode owns the advantage that both J/tp and photon can be tagged cleanly. 
The presence of J / ip is particularly helpful to reduce the combinatorial background. In this 
regard, this mode is practically much more useful than the purely electromagnetic decay 
Vb -> 77- 

By an explicit pQCD calculation based on NRQCD approach, we infer the branching 
ratio of this process to be of order 10~ 7 . Although the absolute value of this ratio is small, it 
is already larger than that of the hadronic decay mode rjb — > J/ipJ/ip, which was previously 
thought of as a golden mode for searching 775. Our analysis indicates that the chance of 
observing rjb through this decay channel, followed by J/ip decay to a lepton pair, with a 
corresponding branching ratio about 10~ 8 , seems still open in Tevatron Run 2, and is quite 
promising in the forthcoming LHC experiment. However, one should bear in mind that 
the major QCD background events are expected to greatly outnumber the desired signal 
events. A thorough study of the associated J/ip + 7 production is welcome, in order to help 
experimentalists to impose optimal kinematical cuts to singlet out the signal events from 
the abundant background events. The transverse polarization of J/ip in the signal events 
should be employed to effectively veto the background. 

While this paper is being written, we are informed of the same rjj, — > J/-07 process 



15 



being also considered by Gao, Zhang and Chao [30j]. These authors evaluated loop integrals 
numerically. It has been checked that once the same input parameters are assumed, our 
numerical prediction for the branching ratio is compatible with theirs. These authors have 
also studied various other radiative decay channels of bottomonium to charmonium. 
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APPENDIX A: DERIVING ANALYTICAL EXPRESSION FOR / 



In this Appendix we illustrate how to reduce the one-loop four-point function in ([5} to the 
sum of much simpler two- and three-point scalar integrals. We start with the dimensionless 
integral /: 

(k-Q)q 2 -(k-q) (Q-q) 



f 



\m 2 ) 



d*q- 



(Al) 



(q + Q) 2 (q - Q) 2 (q 2 + 2k-q- P 2 ) (q 2 - 2k ■ q - P 2 ) ' 

where Q 2 = 4m| and P 2 = Am 2 . The +ie prescription in the propagators is implicitly 
implied. 

First note the second term in the integrand can be simplified by using the identity of 
fractional sum: 



d A q 



(k -q)(Q- q) 



(q + Q) 2 (q - Q) 2 (q 2 + 2k- q - p 2 ) ( q 2 -2k- q - p 2 ) 



16 ' ' 



l(q-Q) 2 (q + Q) 2 



1 



(q ~ Q) 2 



q 2 - 2k ■ q - P 2 q 2 + 2k-q-P 2 
1 



(A2) 



_q 2 ~2k-q-P 2 q 2 + 2k-q-P 2 _ 

which is nothing but the scalar 2-point functions, thus can be trivially worked out. 

Disentangling the first term in (1A1I) needs slightly more labor. Since the integrand is an 
even function of q, we are free to add terms linear in q in the numerator, without influencing 
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the result: 



(q + Q? (q - Q) 2 (q 2 + 2k-q- P 2 ) (q 2 -2k-q- P 2 ) 

I /\ 4 {Q 2 + P 2 )q 2 + 2Q 2 k-q + 2P 2 Q-q + Q 2 P 2 -Q 2 P 2 
a q 



+ 



Q 2 + P 2 

Q 2 
Q 2 + P 2 
P 2 



d"q 
d 4 q 



(q + Q) 2 (q - Q) 2 (q 2 + 2k- q - p 2 ) ( q 2 -2k- q - p 2 ) 
1 



(q + Q) 2 (q - Q) 2 (q 2 - 2k ■ q - p 2 ) 

1 



(A3) 



Q 2 + P 2 J * (q - Q) 2 (q 2 + 2k-q- P 2 )(q 2 - 2k ■ q - P 2 ) ' 

which consists of two independent 3-point scalar integrals. They can also be worked out in 
closed form, following the method outlined in Ref. 29]. 

With the aid of ( 1A2I) and ( 1A3I) . we can decompose the original / into three pieces: 



f{u) = /1(U)+/2(U)+/ 3 (U), 



(A4) 



where u = P /Q = mym 2 , and 



7T- 1 



d 4 q 



(q-Q) 2 



1 



In 



/ a («) 



2 (1 - u) 

2-u 

8Q 2 k-Q 
tn 2 {Q 2 + P 2 ) 



2(1 -u 
d^q 



q 2 -2k-q-P 2 q 2 + 2k-q-P 2 



2 — it 



(A5) 



1 



(g + Q) 2 (? - Q) 2 (q 2 - 2k ■ q - p 2 ) 

u 



2 \ I 

< In 2 2 H — In 2 -u + ln]2 - u] In 

1+u I 2 L J 



2(1 -u) 



h{u) 



+ 2Li 2 [-«] + Li 2 

8 P 2 k ■ Q 

~ m 2 {Q 2 + P 2 ) 

2u , , 
in 



it — 1 
2u 

d 4 q 





~u~ 




u 2 — u 




+ 2Li 2 




- Li 2 




} 




.2. 


2 





+ In u In 

2yri 



1 + u 



l-u_ 
ln[2 - u] , (A6) 



(q ~ Q) 2 (q 2 + 2k-q- P 2 )(q 2 - 2k ■ q - P 2 ) 



1 + u 
+ 2ni 



u 



it 



1 + U 



2-u 
In 



In 



a 



u 



2(1 -u) 



+ Li 2 



2 - 



2-u 



(A7) 



One can then readily reproduce the analytic results shown in (jHJ) and (ITU]) . 



[1] For a comprehensive review on various aspects of see N. Brambilla et al, CERN-2005-005 
(CERN, Geneva, 2005) [arXiv:hep-ph/0412158| , and references therein. 



17 



[2] D. Ebert, R. N. Faustov and V. O. Galkin, Phys. Rev. D 67, 014027 (2003) 
[arXiv:hep-ph/021038l] . 



S. Recksiegel and Y. Sumino, Phys. Lett. B 578, 369 (2004) |arXiv:hep-ph70 305178]. 

B. A. Kniehl, A. A. Penin, A. Pineda, V. A. Smirnov and M. Steinhauser, Phys. Rev. Lett. 



92, 242001 (2004) [arXiv:hep-ph/ 0312086|. 

A. Gray, I. Allison, C. T. H. Davies, E. Gulez, G. P. Lepage, J. Shigemitsu and M. Wingate, 
Phys. Rev. D 72, 094507 (2005) |arXiv:hep-lat/0507013| . 

A. Heister et al. [ALEPH Collaboration], Phys. Lett. B 530, 56 (2002) [arXiv:hep-ex/0202011] . 

M. Levtchenko [L3 Collaboration], Nucl. Phys. Proc. Suppl. 126, 260 (2004). 

J. Abdallah [DELPHI Collaboration], Phys. Lett. B 634, 340 (2006) |arXiv:hep-ex/0601042| . 

M. Artuso et al. [CLEO Collaboration], Phys. Rev. Lett. 94, 032001 (2005) 

[arXiv:hep-ex/0411068| . 

E. Braaten, S. Fleming and A. K. Leibovich, Phys. Rev. D 63, 094006 (2001) 
|arXiv:hep-ph/0008091| . 

J. Tseng [CDF collaboration], FERMILAB-CONF-02-348-E, Presented at 5th International 
Conference on Quark Confinement and the Hadron Spectrum, Gargnano, Brescia, Italy, 10-14 
Sep 2002. 

F. Maltoni and A. D. Polosa, Phys. Rev. D 70, 054014 (2004) [arXiv:hep-ph/0405082| . 



Y. Jia, arXiv:hep-ph/0611130 



D. P. Roy and K. Sridhar, Phys. Lett. B 341, 413 (1995) |arXiv:hep-ph/ 9407390 1 . 

C. S. Kim, J. Lee and H. S. Song, Phys. Rev. D 55, 5429 (1997) |arXiv:hep-ph/9610294| 



P. Mathews, K. Sridhar and R. Basu, Phys. Rev. D 60, 014009 (1999) (arXiv:hep-p h/9901276 |. 
G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev. D 51, 1125 (1995) [Erratum-ibid. D 



55, 5853 (1997)] |arXiv:hep-p h/9407339|. 

N. Brambilla, D. Eiras, A. Pineda, J. Soto and A. Vairo, Phys. Rev. D 67, 034018 (2003) 
[arXiv:hep-ph/0208019] . 

N. Brambilla, Y. Jia and A. Vairo, Phys. Rev. D 73, 054005 (2006) | |arXiv:hep-ph/0512369| . 
S. Fleming, C. Lee and A. K. Leibovich, Phys. Rev. D 71, 074002 (2005) 
[arXiv:hep-ph/0411180] . 

B. Guberina and J. H. Kuhn, Lett. Nuovo Cim. 32, 295 (1981). 

J. G. Korner, J. H. Kuhn, M. Krammer and H. Schneider, Nucl. Phys. B 229, 115 (1983). 



18 



[23] W. M. Yao et al. [Particle Data Group], J. Phys. G 33, 1 (2006). 

[24] V. V. Braguta, A. K. Likhoded and A. V. Luchinsky, Phys. Rev. D 72, 094018 (2005) 
|arXiv:hep-ph/0506009] . 

[25] J. P. Ma, Nucl. Phys. B 605, 625 (2001) [Erratum-ibid. B 611, 523 (2001)] 

[arXiv:hep-ph/0103237] . 
[26] Y. D. Yang, |a?XT^:hep-ph/ 0404018 

[27] G. Li, T. Li, X. Q. Li, W. G. Ma and S. M. Zhao, Nucl. Phys. B 727, 301 (2005) 
[arXiv:hep-ph/0505158] . 

[28] Y. J. Gao, Y. J. Zhang and K. T. Chao, Chin. Phys. Lett. 23, 2376 (2006) 

|arXiv:hep-ph/0607278] . 
[29] G. 't Hooft and M. J. G. Veltman, Nucl. Phys. B 153, 365 (1979). 
[30] Y. J. Gao, Y. J. Zhang and K. T. Chao, hep-ph/0701009[ 



19 



